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At the microscopic level, equilibrium liquid's translational symmetry is spontaneously broken at 
the so-called dynamic glass transition predicted by the mean-field replica approach. We show that 
this fact implies the emergence of Goldstone modes and long-range density correlations. We derive 
and evaluate a new statistical mechanical expression for the glass shear modulus. 
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Introduction. — The difference between a liquid and 
a solid manifests itself in their response to a volume pre- 
serving shear deformation: a liquid flows whereas in a 
solid internal forces arise which are proportional to and 
oppose the deformation. Solid's rigidity is an emergent 
[l| property: although it amounts to transmitting inter- 
nal forces opposing a volume preserving deformation over 
macroscopic distances, it does not originate from long- 
range intermolecular interactions. However, it is intu- 
itively plausible, and it was argued at a phenomenolog- 
ical level [3, 01 , that the emergence of rigidity has to be 
accompanied by appearance of long-range correlations. 

Interestingly, the first microscopic argument for the ex- 
istence of long-range correlations in crystalline solids was 
given in 1966 by Wagner (3], preceding phenomenological 
analyzes of Refs. [ITSl- However, the part of his paper 
dealing with solids seems to have been overlooked. More 
importantly, while he showed that in crystalline solids 
broken translational symmetry implies long-range corre- 
lations, he did not connect the latter to solids' rigidity. 

The intimate relation between long-range correla- 
tions and rigidity was established somewhat indirectly. 
Bavaud et al. [5[ analyzed the standard expression for 
the elastic moduli, which is derived by calculating free 
energy change due to a deformation. They proved that 
this expression reduces itself to the inverse of the isother- 
mal compressibility if both intermolecular interactions 
and density correlation functions are short range. Thus, 
they showed that for systems with short-range interac- 
tions the absence of long-range density correlations im- 
plies the vanishing of the shear modulus. 

The microscopic argument for the existence of long- 
range correlations in crystals was independently derived 
by Szamel and Ernst [6|. They proposed an expression 
for a displacement field in crystals in terms of the mi- 
croscopic density and showed that the displacement held 
defined in this way exhibits long-range correlations, as ex- 
pected on the basis of phenomenological arguments 

Here we consider an amorphous solid which appears as 
a consequence of the so-called dynamic glass transition 
predicted by the mean-field replica approach We 
show that, in close analogy with crystalline solids, broken 
translational symmetry of this amorphous solid implies 
the emergence of long-range correlations and rigidity. 



It should be recalled that, although macroscopically 
the amorphous solid's density is constant, microscopi- 
cally it is randomly non-uniform. Thus, a rigid transla- 
tion of the amorphous solid produces an equivalent but 
different state. To describe this randomly broken trans- 
lational symmetry, and amorphous solids ingeneral, we 
use the replica approach of Franz and Parisi [30 . In this 
approach, the appearance of the amorphous solid man- 
ifests itself in non-trivial replica off-diagonal densities. 
We argue that the existence of distinct, rigidly shifted 
states of the amorphous solid manifests itself in a family 
of equivalent sets of replica off-diagonal densities. Next, 
we propose an expression for the displacement field in 
terms of microscopic two-point replica off-diagonal den- 
sities. We show that the displacement field defined in this 
way exhibits long-range correlations. Since the displace- 
ment field is defined in terms of microscopic two-point 
replica off-diagonal densities, long range of the displace- 
ment field's correlations implies a slow decay of a four- 
point replica off-diagonal correlation function. Finally, 
we derive and evaluate a new expression for the shear 
modulus of amorphous solids. To this end we evaluate 
at the microscopic level the force needed to maintain a 
shear deformation and compare it to an expression ob- 
tained from the macroscopic theory of elasticity. 

Family of equivalent glassy states. — Following Franz 
and Parisi [8|, |9| we consider a system of particles coupled 
to a frozen (referred to as "quenched") configuration of 
the same system via an attractive potential 
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where and are positions of particles i and j in the 
system and in the quenched configuration, respectively, 
and w{r) is a monotonic function with a minimum at 
r — 0. We assume that when the amplitude e of the 
potential is gradually reduced, if the temperature is low 
enough or if the density is high enough, the system gets 
stuck in a metastable state characterized by non-trivial 
correlations with the quenched configuration and that 
this metastable state survives in the e limit. The 
appearance of such a metastable state is termed the dy- 
namic glass transition To facilitate averaging over 
a distribution of quenched configurations we replicate the 
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system s times. Including the quenched configuration as 
the 0th replica we have m = s+1 times replicated the sys - 
tem and at the end we take the limit s — ?■ or to — >■ 1 . 
At the level of the replicated system, the dynamic glass 
transition manifests itself in the appearance of non-trivial 
two-point replica off-diagonal densities nap{r), a ^ fi. 
With w(r) having a minimum at r = 0, these functions 
exhibit a pronounced peak at r — 0, oscillate in phase 
with the equilibrium pair distribution function g(r) and 
decay to r? at large distances (n is the average density) . 
In the following the metastable state with all replica off- 
diagonal densities having the main peak at r = will be 
referred to as the "classical" state. 

The above described construction can also be per- 
formed with a potential that pins the system in a position 
shifted with respect to the quenched configuration, 
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The metastable states obtained with potentials ([T]) and 
© are identical up to a rigid shift by vector a of the 
system relative to the quenched configuration. At the 
level of the replicated system this translation corresponds 
to the following transformation of two replica densities, 

»^ao(?') -> "Qodr - a|), noa('') ^ ^loadr + a|) (3) 

for a > and keeping all other two replica densities (both 
diagonal and off-diagonal) unchanged. In the following 
we consider only densities n^o, a > 0, since densities riQa 
can be obtained from the former ones. 

This observation implies that instead of one glassy 
state we have a continuous family of states that can be 
labeled by a vector a denoting translation with respect to 
the "classical" state. This is analogous to what is found 
in other systems with broken continuous symmetry, e.g. 
in ferromagnets Q or in crystals 

We note that the rigid shift described above breaks 
replica symmetry. However, since replica non-symmetric 
states obtained by rigid translations arise due to broken 
translational symmetry, they have the same free energy 
as the "classical" solution (which is replica symmetric 

0). Thus, the present case is fundamentally different 
from replica symmetry breaking in mean-field spin glasses 
(l^ . In the latter case, replica symmetry breaking origi- 
nates from the existence of a multitude of different glassy 
states unrelated by translations. Importantly, the replica 
non-symmetric state minimizes the free energy. 

Microscopic definition of displacement field. — Ac- 
cording to the conventional definition, a displacement 
field in crystalline solids is defined in terms of depar- 
tures of particles' instantaneous positions from their lat- 
tice sites. Notably, this definition implicitly assumes that 
each particle can be assigned to a specific lattice site. 
Thus, it is not applicable for, e.g. crystals with vacancies 
or interstitials. This fact was one of the motivations for 
a microscopic definition of the displacement field [l^ . 



Since amorphous solids do not have any underlying 
crystalline lattice, it is impossible to use the conventional 
definition. An operational definition that is closest to the 
one used for crystals is as follows. First, average positions 
of particles are identified. Then, the displacement field 
is defined in terms of departures of the instantaneous 
positions with respect to the average positions. While 
a procedure like this can be implemented in a computer 
simulation or an experiment 15 1, it is not clear how 
to formulate it theoretically. The definition we propose 
is similar in spirit to this procedure in that it relates 
microscopic and average densities of the amorphous solid. 

We propose the following definition of (the Fourier 
transform of) the displacement field in terms of micro- 
scopic replica off-diagonal density. 
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Here nQo(ri,r2) is the "classical" two-point replica off- 
diagonal density and J2i j ~ J*" )^(r2 — r°) is the mi- 
croscopic two-point replica off-diagonal density with rf 
and being positions of particles i and j in replica 
a and 0, respectively. Finally, Af is the normalization 
factor, M = (3s)~i / dr2i J2a>a (Sri'^aolri, ra))^. Note 
that definition Q is symmetric with respect to replica 
indices a, a > 0. This reflects the fact that the system is 
first deformed and only subsequently replicated. Thus, 
deformations in all a > replicas are the same. 

Importantly, definition (j4]) is only applicable to in- 
finitesimally small deformations of the "classical" state. 
A similar restriction applies to the definition of the dis- 
placement field in crystals proposed in Refs. p, |l3|. 

Long-range density correlations. — First, we adapt 
the calculation presented in Ref. ^ and prove that cor- 
relations of the displacement field defined through Eq. 
(m exhibits a small wavevector divergence. Then we 
show that long-range displacement field correlations im- 
ply slow decay of a component of a four-point replica 
off-diagonal correlation function. 

We start with Bogoliubov's inequality 0, 

> I {AB) |2 with A = F-i/2n . u*(k) and 

B = V^^^^s^^J2a>o^ ' i^aSa(k.) whcrc n is an arbi- 
trary unit vector, go, is the Fourier transform of the 
microscopic momentum density in replica a, ga(k) = 



TOV"e 



, with TO denoting the particle's mass and 



vf being the velocity of particle i in replica, and is 
the Liouville operator Q in replica a. 

The cross term can be evaluated using the self-adjoint 
property of In the k — ;> limit we obtain (AB) = 
—kBT/s. can be expressed in terms of the corre- 

lation function of the stress tensor. The stress tensor in 
replica a is defined through the continuity equation for 
the momentum in replica a, iCagaO^) — CTq, (k;t). 
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Substituting (AB) and (|-Bp) into Bogoliubov's in- 



where k — k/fc. 

Bound (O for the small wavevector behavior of the 
displacement field correlations is compatible with the re- 
sult of a phenomenological analysis which predicts a 
divergence i^j- This agreement, and the closeness of the 
steps outlined above and the general scheme presented in 
Chapter 7 of Ref. ^ provide an additional justification 
for definition (|4|) and for recognizing the displacement 
field as a "hydrodynamic Goldstone mode" . 

Finally, since the displacement field is defined in terms 
of microscopic two-point replica off-diagonal densities, 
the left-hand-side (LHS) of Eq. ([5]) can be re-written 
in terms of four-point replica off-diagonal densities, 

TTTQ 1 /"j J - anio(ri,r2), c'nio(r3,r4) 

X (nioio(ri,r2,r3,r4) - nio2o(ri, r2, rs, r4)) e"*'^''"i^(6) 

The combination of Eqs. ([5]) and ^ proves the slow 
decay of a component of a four-point replica off-diagonal 
correlation function of the amorphous solid. An anal- 
ogous calculation for a crystal proves the slow decay 
of correlations of high Fourier components of the (non- 
translationally invariant) two-particle density 0, 0, Il3j . 

Shear modulus. — The derivation of a new formula 
for the amorphous solid's shear modulus proceeds in 
two steps. First, we consider the change of replica off- 
diagonal densities upon an infinitesimally small, long 
wavelength deformation. Next, we calculate microscop- 
ically the force needed to maintain such a deformation. 
Comparison of this force with the corresponding macro- 
scopic expression allows us to identify the shear modulus. 

According to Eq. ([3]), under an infinitesimally small 
uniform translation, replica ofl[-diagonal densities Uao, 
a > 0, change in the following way: 

naa{r) riaoir) - a- drUaoir). (7) 

To generalize Eq. ([7|) to an infinitesimally small non- 
uniform translation we imagine a pinning potential which 
deforms the system by imposing an infinitesimal shift 
that depends on the position of the particle. This sug- 
gests the following generalization of Eq. ([7]) , 

nao{ri2) nao(i"i,i"2) 

= naQ{ri2) - a(ri) • drj^nao{ri2), (8) 

where ri2 = |ri — and a(ri) is an infinitesimal, long- 
wavelength deformation imposed on the system. We em- 
phasize that ([5]) is only assumed for a long wavelength 



equality, and taking the s — ;> limit we obtain 
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volume-preserving deformation. Specifically, we assume 
a transverse deformation such that dr ■ a(r) = 0. 

Eq. ([U is inspired by the assumption made by Triezen- 
berg and Zwanzig 17], and by Lovett et al. [18] in their 
analyzes of surface tension. A similar assumption was 
made by Szamel and Ernst @ in their derivation of shear 
modulus of crystalline solids. Finally, assumption ([S]) is 
consistent with our microscopic definition of the displace- 
ment field: if all densities riao, a > 0, are changed as in- 
dicated in ([8|) then the Fourier transform of the average 
displacement field is equal to (u(k)) — J dre^'^^''^a{r). 

Next, we calculate the force exerted on one replica of 
the system by the pinning potential chosen in such a way 
that it imposes the replica off-diagonal densities ([SJ. 

First, we identify the pinning potential needed to main- 
tain densities Since we are only interested in in- 
finitesimally small deformations and thus infinitesimally 
small changes of replica off-diagonal densities, we resort 
to a (static) linear-response type relation. Thus, we ex- 
press the pinning potential using a functional derivative 
of the inter-replica potential with respect to the replica 
off-diagonal density, 

V / drsdri ( i^""'f^'^^l ) [-a(r3) • dr^Ufjoir^i)] . 
f^J V^"^o(r3,r4)y„ 

(9) 

Note that since we have assumed a volume-preserving 
deformation, the functional derivative in expression Q 
should be taken at constant density. 

Then, we calculate the force per unit volume exerted 
by the above pinning potential on the replica a. The 
Fourier transform of this force, after integrating by parts 
and using the transverse character of the deformation, 
can be written in the following way: 

Fa(k) = / rfri...dr4e-^'' '-- {d,,n^o{ri2)) (10) 

V'^"/3o(r3,r4)/„ 

Finally, we take the long wavelength limit and expand 
the right- hand- side of Eq. ([10]) in a power series in k. The 
zeroth order term vanishes because no force is needed in 
order to accomplish a rigid shift of the amorphous solid. 
The first order term vanishes by symmetry. The second 
order term, in the s — limit and for a shear deformation 
involving y-dependent translation along the x axis, can 
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be used to evaluate the shear modulus of the hard sphere 
glass whereas the standard formula is not applicable for 
hard sphere systems. In contrast to the result obtained in 
a recent replica approach study [llj (which was based on 
the standard formula) we found a discontinuous change 
of the shear modulus at the dynamic glass transition. 
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FIG. 1: Dimensionless shear modulus, fia"^ /ksT, for the hard 



sphere glass as a function of the volume fraction. 
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The shear modulus changes discontinuously at the dynamic 
transition which occurs at (he — 0.5159. 



be written in the following form 
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where /i is given by 
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J dri...dr4 yj^ {dxiUioiru)) {dx^nioirsi)) 



(5nio(r3,r4) V <5n2o(r3,r4) 



Macroscopically, for an isotropic solid the force needed 
to maintain a long wavelength shear deformation is given 
by a formula identical to Eq. (Ilip . where fj, is the shear 
modulus. This fact allows us to identify fi given by Eq. 
(IT2]) as the shear modulus of the amorphous solid. 

To evaluate fi one has to use a concrete (albeit nec- 
essarily approximate) implementation of the replica ap- 
proach in order to calculate replica off-diagonal densities 
and functional derivatives in Eq. (|12p . To obtain Uaoi''') 
we use the recently proposed version [lij which is con- 
sistent with mode-coupling theory [2^. To calculate the 
derivatives we use an additional approximation. 



n'^Scao{ri,r2) = -nQo(i"i, i"2)/3VQo(ri, r2)- 



Here 6cao is the change of the replica off-diagonal direct 
correlation function due to a weak replica off-diagonal po- 
tential Vao- Approximation (jl3p combined with replica 
Ornstein-Zernicke equations allow us to get explicit ex- 
pressions for functional derivatives that enter Eq. (fT2|). 

In Fig. [1] we show the shear modulus for the hard 
sphere glass calculated using as an equilibrium input the 
Percus-Yevick structure factor. 

Discussion. — We identified Goldstone modes and 
long-range correlations appearing in the amorphous solid 
due to the spontaneously broken translational symmetry. 
We derived a new expression for the shear modulus of 
the amorphous solid. This expression is complementary 
to the standard one In particular, our expression can 
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